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1.  INTRODUCTION 

1.1.  Background  and  Overview 

The  theory  and  associated  software  have  previously  been  developed  to  compute  the 
reflectivity  and  transmissivity  of  planar  periodic  structures,  including  active  reflection 
coefficient  of  phased  array  antermas,  using  the  hybrid  finite  element  method  (HFEM)  with 
periodic  boundary  conditions  [1].  Among  its  requirements  was  that  the  entire  non-fi'ee-space 
material  structure  of  the  antenna  or  scatterer  be  contained  within  the  finite  element  region, 
with  the  periodic  radiation  condition  being  applied  on  planar  surfaces  on  or  outside  that 
structure.  This  requirement  is  very  undesirable  in  a  number  of  important  circumstances  in 
which  the  structure  includes  uniform  dielectric  layers.  For  example,  consider  the  case  of  a 
waveguide  array  anteima  that  includes  a  dielectric  sheath  outside  the  array  for  purposes  of 
wide-angle  impedance  matching,  or  a  fi-equency  selective  surface  whose  outer  layers  are 
uniform  dielectrics.  It  was  previously  necessary  to  extend  the  finite  element  mesh  into  the 
dielectric  layers,  with  a  consequent  increase  in  the  size  of  the  matrix  to  be  solved. 

This  report  discusses  a  method  for  incorporating  the  effects  of  uniform  external 
dielectric  layers  into  the  periodic  radiation  condition.  The  technique  was  developed 
previously  for  waveguide  phased  array  analysis  by  multimode  methods  [2],  and  has  also  been 
used  for  fi*equency  selective  surface  analysis  by  moment  methods  [3].  It  is  used  to  modify 
the  matrix  entries  associated  with  the  outermost  surface  of  the  non-uniform  material  structure 
by  finding  an  effective  admittance  that  replaces  the  modal  admittance  of  the  Floquet  modes 
used  to  represent  the  exterior  solution.  This  provides  three  important  advantages:  (1)  the 
size  of  the  matrix  is  reduced  significantly,  especially  when  the  exterior  layers  are  thick  and/or 
dense;  (2)  very  thin  layers  are  permitted;  and  (3)  it  extends  the  range  of  problems  that  may 
be  solved  without  using  periodic  unit-cell  side  wall  boundary  conditions  that  place  demands 
on  the  mesh  that  some  CAD  software  packages  cannot  satisfy.  These  advantages  appear  to 
come  at  very  little  cost  since  the  only  penalty  is  a  few  additional  algebraic  steps  during  the 
matrix  fill  computations. 
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Section  2  summarizes  the  original  theory  without  exterior  dielectrics  and  the 
necessary  modifications.  Section  3  presents  calculated  results  with  comparisons  to  published 
data  obtained  by  other  methods.  Those  results  verify  that  the  additions  to  the  finite  element 
solution  are  effective  for  both  array  antennas  and  fi-equency  selective  surfaces  (FSSs). 

1.2.  General  Problem  Description 

Figure  1  illustrates  the  two  important  classes  of  problems  that  this  theory  addresses. 
Fig.  la  shows  a  side  view  of  an  array  antenna,  each  of  whose  radiating  elements  is  fed  by  a 
single  waveguide.  The  power  incident  fi-om  each  waveguide  is  the  same,  but  the  phase  may 
vary  progressively  fi'om  element  to  element  to  electronically  steer  the  antenna  beam.  The 
answers  that  an  array  designer  may  require  are  the  input  impedance  of  an  element  in  the  array 
environment,  called  the  "active  impedance,"  the  gain  per  element,  or  "active  element  gain," 
and,  if  grating  lobes  are  present,  the  angle  and  power  in  each  lobe.  Accurate  calculation  of 
these  quantities  requires  a  rigorous  solution  that  accounts  for  the  mutual  coupling  between 
array  elements.  Efficient  calculation  relies  on  the  assumption  that  the  array  is  large,  so  that 
edge  effects  are  negligible.  When  that  is  true,  the  analysis  may  be  restricted  to  a  single  array 
unit  cell. 

Fig.  la  shows  that  there  may  be  uniform  dielectric  layers,  or  "superstrates"  above  the 
non-homogeneous  structure.  The  radiation  boundary,  denoted  Fb  may  be  placed  anywhere 
above  that  part  of  the  structure  that  is  non-uniform.  As  indicated  in  the  figure,  Fg  may  cut 
through  a  dielectric,  as  long  as  that  part  of  the  dielectric  above  Fg  is  a  continuous,  planar 
layer. 

Figure  lb  shows  a  side  view  of  a  generic  planar  periodic  structure.  The  excitation  is 
a  plane  wave  impinging  on  the  lower  side.  The  designer  may  need  to  know  the  amount  of 
power  reflected  and  transmitted,  both  co-polarized  and  cross-polarized.  Also,  if  there  are 
grating  or  Bragg  lobes,  their  directions  and  relative  power  and  phase  may  be  important.  As 
the  figure  indicates,  there  may  be  continuous  dielectric  layers  both  above  and  below  the 
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Figure  1.  Generic  Problem  Classes:  (a)  Phased  Array  Antenna  with 
Waveguide  Feeds;  and  (b)  Planar  Periodic  Structure 
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structure  that  are  outside  the  analysis  region,  which  is  formed  by  the  unit  cell  walls  and  the 
upper  and  lower  radiation  boundaries. 

Figure  2  defines  some  of  the  notation  used  throughout  this  report.  The  structures  are 
considered  to  be  infinite  in  two  directions  and  oriented  perpendicular  to  the  z  axis.  The  array 
lattice  may  have  an  arbitrary  skew  angle  y.  The  unit  cell  boundary  is  not  necessarily  a 
parallelogram  as  illustrated,  but  may  have  protrusions  on  one  side  and  corresponding 
indentations  on  the  other,  as  long  as  the  -x  and  +x  side  walls  are  identical  except  for  a 
translation  of  d„  and  tiie  -y  and  +y  side  walls  are  identical  except  for  a  translation  of  dyin  the 
y  direction  and  dyC0X{Y)  in  the  x  direction. 

y 


Figure  2.  Array  Lattice  Geometry  and  Notation 
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2.  THEORY 


2.1.  The  Weak  Form  of  the  Wave  Equation 

The  finite  element  solution  is  based  on  discretization  of  the  “weak  form”  of  the  wave 
equation.  The  wave  equation  for  the  electric  field  in  linear,  isotropic,  inhomogeneous  media 
is 


Vx  Vx£  £  =  0 


(1) 


where  ko  is  the  free  space  wavenumber,  and  the  relative  permeability,  and  permittivity, 
e„  may  be  functions  of  position.  The  inner  product  of  (1)  with  a  weighting  function  W  is 


f  w*' 

Vx  —  VxE  -  kn  e  E 

J 

Q 


dv  =  a 


(2) 


Using  a  Green's  identity: 


/ 


J-VxJF**Vx£  -ito 


€  W  *'E 

r 


dv 


W  n  ds 


=  0 


(3) 


where  F  is  the  union  of  all  boundaries  enclosing  Q  as  well  as  conducting  surfaces  inside  Q. 
Equation  (3)  is  called  the  "weak  form."  The  wave  equation  is  satisfied  in  an  average,  or 
weak  sense.  Note  that  the  Green's  identity  shifted  one  of  the  derivatives  from  E  to  W, 
weakening  the  differentiability  requirement  on  E,  which  will  permit  approximation  of  the 
field  with  fimctions  that  only  have  one  derivative  (linear).  The  final  form  is  obtained  by 
using  the  Maxwell  curl  equation  VxE^  the  constitutive  relation  cofi  =  kofj: 


/ 

Q  L 


— Vx  W  **Vx£ 


kl  € 

0  r 


W  ‘E 


-jk^  rjg  j  W  (n  xH)  ds  =0 


(4) 
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The  term  on  the  right,  a  surface  integral  over  all  boundaries  enclosing  Q,  provides  the  mecha¬ 
nism  for  enforcing  boundary  conditions.  On  the  open  parts  of  F,  suitable  expressions  for  the 
tangential  magnetic  field  may  be  substituted  to  impose  radiation  boundary  conditions.  The 
main  objective  of  this  report  is  to  document  those  expressions  that  are  suitable  for  exterior 
dielectric  layers.  The  expressions  for  the  interior  finite  element  matrix  arising  from  the  first 
integral  in  (4)  are  adequately  described  in  [1]  and  will  not  be  repeated  here.  However,  the 
earlier  derivation  for  the  boundary  terms  arising  from  the  right  side  integral  in  (4)  will  be 
reviewed  because  they  are  a  necessary  starting  point  for  the  present  extension  to  the  theory. 


2.2.  Mode  Functions  and  the  Periodic  Integral  Equation 


The  top  and  bottom  of  the  unit  cell  are  free  space  boundaries  where  a  radiation 
condition  must  be  applied.  The  method  for  doing  so  is  to  substitute  an  integral  equation  for 
« 'xR  into  the  boundary  term  of  the  weak  form.  Due  to  the  periodic  nature  of  the  problem, 
the  integral  equation  may  be  written  as  a  summation  over  "Floquet  modes"  [2]. 

Let  the  electric  field  above  and  below  the  structure  be  expanded  in  a  set  of  mode  func¬ 
tions  gp„,„  where  the  subscript  is  1  for  transverse  electric  (i.e.  E±z)  or  2  for  transverse 
magnetic  (i.e.  H^z): 


g 


pmn 


^  •'  '■  xntfi  ' 

pmn 


(5) 


(xk  -yk  ) 

y  xmn  ymn 

(f  sin  -y  cos(^^) 


w ^  0 ,0 

m,n  =0,0 


(6) 
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2mn 


Jd  d 

y  X  y 


(x  k  +yk  ] 

\  xmn  ^  ymn  / 


+k^ 

xmn  ymn 


m,n  7^  0,0 


cos  +j?  sin^^  j  w,«=0,0 


k  =  kf,  sin  d„  cos  (b. 

xmn  0  0  ^0 


2tc/m 


,  1  •  a  •  A.  27i:«  iTZm  cot  y 

^yn,n  =  Slti  0^  SlTi  (p^  -  —  +  - - - L 


(7) 

(8) 

(9) 


where  (^o»^o)  is  the  main  beam  angle,  (m,n)  is  the  mode  index  and  f  is  the  lattice  skew  angle. 
The  index p  is  1  for  modes  whose  electric  field  is  transverse  to  z  (TE),  or  2  for  modes  whose 
magnetic  field  is  transverse  to  z  (TM).  Note  that  the  mode  functions  are  independent  of  the 
medium’s  constitutive  properties-the  same  functions  will  be  used  to  represent  fields  in  all 
dielectric  media  as  well  as  those  in  fi-ee  space. 

At  any  point  above  the  structure,  the  transverse  electric  field  may  be  expanded  in 
modes  traveling  in  the  +z  direction: 

E'(z^O)  =  i:  E  E  c;„„  (10) 


K 


=  ± 


kill  € 

\)  *  r  r  xmn  ymn 


1/2 


(11) 


where  Cp„,„  is  an  unknown  coefficient.  The  sign  of  the  propagation  constant  k„„  is  chosen 
as  positive  for  outgoing  waves  and  negative  for  incoming  waves  [2:43].  The  transverse 
magnetic  field  in  each  mode  is 
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fx/f'  (z>0)  =  -Y  e' 

pmn  ^  '  pmn  pm 


(12) 


where  Yp„„  is  a  modal  admittance: 


Y 


\mn 


(13) 


Imn 


k.  € 

0  r 
mn  *0 


(14) 


The  total  magnetic  field  above  the  structure  is 

fx^'(z>0)  =  -  E  E  E  ^ImnY (15) 

pmn 


In  order  to  write  the  right  side  in  terms  of  E‘,  the  orthogonality  property 

f  S  ‘Snii  dx  dy  -  6  3.3 

I  ^pmn  ^qij  pq  mt  nj 

unit  cell 

is  used  to  solve  (10)  for  the  coefficients: 

C  *  =  f  E  *  ’  s*  dx  dv 

^  pmn  I  ^  opmn 

unit  cell 


(16) 


(17) 


with  the  result 


fx/f 


'(z^O)  =  -EEE 


^ pmn  ^ pm 


■  / 


unit  cell 


E‘-g* 


pmn 


dx  dy  e 


(18) 


On  the  lower  unit  cell  boundary,  the  transverse  electric  field  includes  an  incident  field 
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term  in  one  of  the  two  lowest  order  modes  depending  on  whether  it  is  transverse  electric 
{q=\)  or  transverse  magnetic  {q=2)\ 

2  oo  00 

^ ^ 0 )=  f,„. E  E  E  c. ('») 

P~\  /«  =  -<»  77  = -oo 

The  signs  of  the  exponential  terms  indicate  the  direction  of  propagation.  Again,  using  the 
modes'  orthogonality  property, 

/  E*  ‘g*„  dx  dy  -  d  S 

pmn  J  ^pmn  pq  mO  nO  (20) 

unit  cell 


giving  the  magnetic  field  below  the  structure  as 


£xH\z<0)  =  ^  ^pmnSpmn  f  ^  ^  '  S  pmn  dy  C  ^ 

pmn  •' 

unit  cell  (21) 

-  2  7  F  e  ^ 

■*  9OO  *g00 


(This  is  a  correction  to  eq.  37  of  [1:13],  which  had  the  signs  reversed.)  The  above 
expressions  (18)  and  (21)  are  the  integral  equations  that  are  substituted  into  the  weak  form 
(4)  to  enforce  the  radiation  conditions.  The  following  section  derives  expressions  for  the 
resulting  matrix  terms. 


2.3.  Periodic  Radiation  Boundary  Terms 

The  matrix  entries  pertaining  mesh  edges  s  and  t  in  the  upper  radiation  boundary, 
are  found  by  substituting  an  ejqjansion  fiinction,  i^(the  “vector  edge  elemenf”),for  in  (18), 
then  substituting  (18)  for  h  xR  and  IJr^  for  W*  in  Eq.  (4): 


Kfs‘S*pmn^dxdye 

*l  n  rn  d 


pmn 


(22) 
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Each  integral  term  contains  a  term  in  exp{±j(A;j„^  +  ky^,y)} ,  which,  besides  the  finite  ele¬ 
ments,  is  the  only  term  with  x  or  y  dependence.  Let  J,(k^,ky)  denote  the  two-dimensional 
Fourier  transform  of  that  part  of  lying  on  the  radiation  boundary: 


/  f  if, 


e  ax  ay 


=  h  •  £{k  ,k  ) 

pmn  xmn^  ymn' 


(23) 


pmn  tm  n 


Now  the  matrix  entries  may  be  written  in  the  more  compact  form 

‘[Y  h  h  ]-^*  (24) 

St  J  ^  *0  ^ tmn  ^  pmn  pmn  pmn  ^  ^ smn  v 

pmn 

Previously,  an  even  more  compact  expression  was  obtained  by  summing  TE  and  TM  modes. 
However,  when  external  dielectrics  are  present,  it  is  important  to  keep  those  terms  separate. 

The  expression  for  the  matrix  pertaining  to  the  lower  radiation  boundary  is  identical 
to  (24).  This  follows  fi’om  the  fact  that  the  summations  in  (18)  and  (21)  are  the  same  except 
for  the  sign  of  the  exponent  inside  the  integral.  In  evaluating  those  integrals,  the  value  of  z 
in  the  exponent  is  taken  as  zero  because  it  represents  only  a  constant  phase  factor. 

When  the  medium  above  and  below  the  finite  element  region  Q  is  a  homogenoeous 
dielectric  half  space  (usually  fi-ee  space),  (24)  may  be  used  directly.  However,  if  the 
dielectric  has  finite  thickness,  or  if  there  are  multiple  dielectric  layers,  then  (24)  must  be 
modified  to  incorporate  an  effective  modal  admittance,  as  described  in  the  following  two 
sections. 
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2.4.  Effective  Modal  Admittance  of  One  Layer 

Consider  first  a  case  in  which  the  +z  radiation  boundary  (Fg)  is  covered  by  a  single 
dielectric  sheet.  Let  7'denote  the  effective  admittance  of  any  Floquet  mode  propagating  in 
the  +z  direction  through  the  sheet,  whose  admittance  for  that  mode  is  Y^.  (The  subscripts 
for  mode  indices  are  dropped  in  the  following  derivation  since  it  applies  equally  to  all  modes. 
The  superscript  A  denotes  air,  while  M  denotes  material.)  At  distance  t  beyond  the  reference 
point,  the  slab  terminates  at  a  half  space  whose  admittance  is  (free  space).  Figure  3  is  an 
equivalent  circuit  depiction  using  a  transmission  line  analogy.  The  reflection  coefficient  at 
the  z  =  ?  boundary  due  to  the  admittance  mismatch  is 


P  = 


Y  j4  M 


z  =  t 

z  =  0 


(a)  (b) 

Figure  3.  Equivalent  Circuit  Representation  of  External  Dielectric  Layer: 
(a)  Geometry;  (b)  Transmission  Line  Model 


Y'  is  the  ratio  HI  Eat  the  reference  point.  The  electric  and  magnetic  fields  carried  by  the 
mode  are  (normalized  to  the  field  strength  of  the  +z  component  at  z  =  0): 
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(26) 


E  = 


e  -  pe 


H  = 


iK^z 


pe 


-jK'^{z*H) 


(27) 


where  is  the  z-directed  propagation  constant  for  the  mode  from  (11).  The  first  term  in 
each  of  (26)  and  (27)  represents  the  wave  traveling  in  the  +z  direction,  while  the  second  term 
represents  the  reflected  wave  traveling  in  the  -z  direction.  Taking  the  ratio  of  magnetic  to 
electric  field  at  z  =  0  and  multiplying  the  numerator  and  denominator  by  exp(/  ^  t) : 


Y  f  —  Y  ^  — 


JK^t 


pe 


-JK^t 


-  pe 


(28) 


which  is  the  usual  expression  for  an  admittance  discontinuity  transformed  back  a  distance 
t  along  a  uniform  transmission  line  whose  characteristic  admittance  is  [4:10].  A  more 
convenient  expression  is  obtained  using  the  Euler  identities: 


y'  =  cos(/r^O  +;(T^/7-^)sin(A:^0 
(7"/7^)cos(a-^0  +7sin(A:"0 


(29) 


An  important  feature  of  Floquet  modes,  like  any  waveguide  modes  that  can  be 
expressed  in  terms  of  orthogonal  transverse  mode  functions,  is  that  any  discontinuity 
transverse  to  the  direction  of  propagation  (z  in  this  case)  does  not  cause  cross-coupling 
between  modes.  The  {m,n)  mode  incident  on  a  dielectric  boundary  parallel  to  the  x-y  plane 
will  reflect  energy  only  into  the  {m,n)  mode  in  the  dielectric,  and  transmit  energy  only  in  the 
{m,n)  mode  in  free  space.  Consequently,  the  earlier  expression  for  boundary  matrix  entries, 
Eq.  (24)  may  be  modified  simply  by  replacing  Yp„„  with 


7 

pmn 


jM 

pmn 


i^p1nJ^pmn}<^os(Kr^„0  +ysin(/r^„0 


(30) 
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where  anA  denote  the  modal  admittances  in  uniform  media,  dielectric  and  free 
space,  respectively. 

Some  caution  is  needed  in  calculating  (30)  numerically  because  the  complex  sine  and 
cosine  involve  sinh  and  cosh  functions.  For  large  tn  and  n,  is  a  large  complex  number 
and  may  cause  floating  point  overflow  in  those  function  evaluations.  Appendix  A  describes 
a  stable  approach  to  evaluating  (30)  that  avoids  the  problem. 

2.5.  Effective  Admittance  for  Multiple  Layers 

Consider  next  a  case  in  which  there  are  two  exterior  dielectric  layers,  with  the 
innermost  having  a  thickness  and  characteristic  modal  admittance  Figure  4(a)  shows 
a  transmission  line  analogy  for  this  situation.  The  second  (outermost)  layer  is  identical  to  the 
one  considered  in  the  previous  section.  Noting  that  the  portion  of  the  circuit  to  the  right  of 
the  second  dotted  line  is  identical  to  that  shown  earlier  in  Fig.  3,  the  reduced  model  shown 
in  Fig.  4(b)  is  equivalent.  Now,  noting  the  similarity  of  Fig.  3(b)  to  Fig.  4(b),  it  is  evident 
that  T"is  given  by  Eq.  (30)  with  the  substitutions  and  Y^  ^  Y".  By 

extension,  Eq.  (30)  can  be  used  as  a  recursion  formula  for  any  number  of  layers. 


(a)  (b) 

Figure  4.  Equivalent  Circuit  Representation  for  Two  External  Dielectric  Layers: 
(a)  Transmission  Line  Model;  (b)  Reduced  Model  Using  Results  of  Section  2.4 
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2.6.  Incident  Field  Terms 


2.6. 1.  Free  Space 

The  right  side  vector  in  the  matrix  equation  arises  from  the  second  term  on  the  right 
side  of  (21),  due  to  a  unit  amplitude  incident  electric  field  arriving  from  z  <  0.  Its  form  is  a 
plane  wave,  which  is  the  0,0-order  Floquet  mode.  The  vector  entry  for  edge  t  is: 

=  -JK%  f  (31) 


The  constant  phase  factor  exp{-y  k;„„  z)  in  (21)  is  ignored.  Taking  all  of  the  terms  outside  the 
integral  that  do  not  depend  on  x  and  y: 


=  -VK 


0  %  ^qOO  ^qOO 


y) 


dx  dy 


(32) 


Note  that  the  integral  is  the  two-dimensional  Fourier  transform  of  the  finite  element's 
tangential  component  on  F^,  which  is  denoted  Note  also  the  following  simplification 
from  (8)  and  (9): 


When  the  medium  is  free  space,  the  dominant  mode  admittance  is  then 


Y 


900 


cos 


% 

1 


%  cos 


q  =  1 
q=2 


(34) 


Now  (32)  may  be  rewritten  in  the  form 
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-VK 


cos  (  X  sin  00  -y  cos  0^ )  ’  Coo  ^  =  1 

sec  00 +j;sin  0o)- Coo  ^=2 


(35) 


2.6.2.  With  Exterior  Dielectric  Layer 

For  the  case  in  which  there  is  a  dielectric  layer  covering  the  side  of  a  periodic 
structure  on  which  an  incident  wave  impinges,  the  incident  field  admittance  term  Y^oo  in  (21) 
will  be  different  than  (34).  To  determine  the  correct  form,  modal  expansions  for  the  fields 
in  and  outside  the  dielectric  will  be  matched  at  the  air-dielectric  interface,  then  the  ratio  7’ 
=  (f  xR)  /  E  will  be  evaluated  at  z  =  0  for  each  mode.  In  the  following  analysis,  the  finite 
element  domain  will  be  assumed  to  occupy  part  of  the  region  z  ^0,  while  the  dielectric  slab 
occupies  the  region  -t  ^  z  <0. 

In  the  fi-ee  space  region  z  <  -t,  the  transverse  electric  field  is  again  represented  by  a 
unit-amplitude  wave  in  the  O'th  order  Floquet  mode  traveling  in  the  +z  direction,  and  some 
unknown  number  of  modes  with  unknown  weights  Cp„„,  traveling  in  the  -z  direction; 


E  (z< -t)  =  g 


Yc 


p,m 


pmn  ^ pmn 


(36) 


which  is  the  same  as  (19)  except  that  the  phase  reference  is  now  the  air-dielectric  interface 
(the  limits  on  summations  over  p,  m  and  n  are  the  same  as  in  (19).  The  corresponding 
magnetic  field  is 


£  'X  H  (z  <  -t)  = 


'^gOO^gOO  ^ 


-JK„Az+t) 


*  E 


Y 

pmn 


c  g 

pmn  ^ 


pmn 


(37) 
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In  the  dielectric  region,  it  is  not  adequate  to  represent  the  +z  traveling  wave  as  the  incident 
mode  only  because  it  may  include  other  modes  scattered  from  z  >  0,  then  reflected  by  the 
dielectric  interface  at  z  =  -t.  Hence,  the  transverse  fields  in  the  dielectric  are  comprised  of 
an  unknown  number  of  modes  traveling  in  each  direction: 


E  ‘  (-t<z  <0)  =  ^ 


g  (A 

^pmn^  pntn 


M 


+  B 


j  IC^„Z 


pmn 


(38) 


Pifn^n 


z^H\-t<z<Q)  =  -  y)  Y^„„g  (A 

^  pmn  ^pmn^  pmn 

P^m^n 


-B 


M 


pmn 


)  (39) 


The  phase  reference  for  the  fields  in  the  dielectric  is  at  z  =  0,  while  that  for  the  fields  in  free 
space  given  by  (36)  and  (37)  is  at  z  =  -t.  This  constant  phase  difference  is  absorbed  into  the 
complex  constants  Ap„„  and  Bp^„. 

Equating  (36)  and  (38)  at  z  =  -t,  then  taking  the  product  with  each  mode  function  and 
integrating  over  the  unit  cell  cross  section,  the  following  result  arises  due  to  the  orthogonality 
of  modes: 


pq  mu  nO 


+  c 


pmn 


=  A  e 

pmn 


B  e 

pmn 


(40) 


Equating  components  firom  (37)  and  (39): 


.  M 


-  B 


pmn 


.  M 

e  ) 


(41) 


Combining  (40)  and  (41)  leads  to  an  expression  for  {Ap„„  -  Bp„„)  in  terms  of  (Ap„„  +  Bp„„)  as 
follows: 
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•  .  .A/ 

A  +  B  e  = 

pmn  pmn 


.  M 


.  M  M  (42) 

pq  mO  wO  y  pmn  pmn  ' 

pmn 


pmn 


-7, 


/?m/2 


^"^pmn  ^pmn^  sin{K^^t)  (43) 


+  2  7  S  d  ^ 

pmn  pq  mO  nO 


(4,„„  -  )  ( 7^^„  cos ( K^„  t)  +;•  7^„„  sm{K^„  t)  )  = 

2  ^pmn  ^pq  ^«0  ^«0  "  pmn  ^p«,„  >(  COS(  Ar^„  0  +y  7^^^  sin(  t)  ) 


(44) 


{A  -  )  = 

pmn  pmn  ' 


2  7 


gOO 


cos  ( /c^  0  +y  7^op  Sin(  /) 

-ri  -i? 

pmn  pmn'  \a  \a  \a 

^DmnCos(A:  0 +y  7  sin(K^„t) 

pmn  V  mn  /  pmn  mn  ' 


(45) 


Referring  back  to  (39),  the  transverse  magnetic  field  evaluated  at  z  =  0  is 


Z^H*  I  =  -  ^nmnZ  "5  ) 

•  pww  ^pmn^  pmn  pmn^ 


2=0 


p,m,n 


(46) 


The  integral  equation  for  the  magnetic  field  in  the  region  z  ^  0  is 
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(47) 


zx/f'(z^0)  =  ^'pmnZpmn  /  ^  ‘  '  S*p„n  dy  e  ^ 


p,m,n 


unit  cell 


-if'  s 

qOO^qOO 


Note  the  primes  on  the  admittanees,  denoting  that  these  are  now  effective  admittances.  In 
addition,  the  tilde  over  Y  in  the  second  term  on  the  right  distinguishes  the  incident  field 
admittance.  Evaluating  (47)  at  z  =  0  with  E‘  from  (38): 

Zxii'  \  =  +  X)  (48) 

r-0  pytn^n 


Equating  this  to  (46)  gives: 


M  - 
pm 


Y  z 

~~tn  o 


(A  -  B  )  =  2 

pmn  ^  pmn  pmn  '  ^00 


-  E  I', 


/  - 


pm 


p,m,n 


P»w 


n  ^pmn  ^^pmn  ^ pmn  )  (49) 


and  exploiting  the  modes’  orthogonality  again: 


(A  -  B  ) 

pmn  pmn  ' 


=  2 


Y 

— 6  6 
M  Pq'nO 

^  pmn 


Y 

_  pmn 
^  pmn 


(A  +  B  ) 

pmn  pmn  ' 


(50) 


Comparing  this  to  (45)  it  is  now  evident  that  the  equivalent  admittance  for  the  incident  field 
term  is 


Y  Y^ 

•*  gOO  ^  qOO 


Y„^n  cos(  a:J1  t)  +  jY 


qOO 


qOO 


•  r  A/  . 

sin(  aTqq  r) 


(51) 
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which  is  to  be  substituted  for  Y^oo  in  (32).  The  earlier  comment  regarding  the  sign  of 
from  (1 1)  is  especially  important  here,  because  it  will  result  in  a  sign  change  in  the  second 
denominator  term. 


2.6.3.  Multiple  Dielectric  Layers 

Suppose  now  that  there  is  a  second  dielectric  layer  occupying  the  space  0  <z<t2,  and 
the  boundary  Fa  isat  z  =  t2  .  The  finite  element  domain  is  now  in  the  region  z  >  /  .  The 
fields  in  the  second  dielectric  layer  may  be  expressed  as  in  (38)  and  (39),  except  with  the 
phase  reference  at  z  =  • 


E\Q<z<t2)  =  5^ 


p^m^n 


^  pmn^ 


a 


.M2  , 


pmn 


+  6 


.M2  ,  . 


pmn 


(52) 


z^H\Q<z<t,)  =  -  V  (a  e  (53) 

2  '  A— ^  ^ pmn^  pmn  pmn 


P>m,n 


Matching  fields  at  z  =  0  with  (38)  and  (39)  and  applying  mode  orthogonality  gives 


M2  .M2 

(A  +B  )  =  (a  b 

^  pmn  pmn^  ^  pmn  pmn  ' 


(54) 


.M2  .M2 

(A  -B  )  =  Y^^  (a  -  b 

pmn  V  pmn  pmn-'  pmn  ^  pmn  nmn  / 


pmn 


(55) 


Combining  the  above  two  with  (50); 
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(56) 


,  .  M2 

if’  6  6  6  -Y  (a  e~^ 

^  ^  qOO  pq  ^mO  ^nO  ^pmn^^pmn^  ^ 

=  7  (a  e  ^  -  b 

^ pmn ^  pmn 


pmn 


pmn 


.  M2 

^  J  ^mn  ^2^ 
.  M2 


which  may  be  solved  for  {ap„„  -  bp„„)  in  terms  of  {ap„„  +  bp„„)  with  the  result: 


(a  -  b  )  = 

pmn  pf^n  ' 


2  y' 

^  ^  ?oo 


M2  .M2  .  .  ^.  /  .,  M2  . 

7^ooCos(;Coo  t^)  sm(/Coo  t^) 

■v'  /  -.  \  ,  •  TT  A/2  .,  M2  . 

-  (a  +  b  )  ^2^  J  ^ pmn  ^2^ 

'^^pmn  ^pmn’  .  M2  .  .^1  M2  ^ 

^ pmn^^^i^mnh^  ^  J  ^ pmn 


(57) 


Note  the  distinction  between  F'  in  the  numerator  and  Y'  in  the  denominator.  The  former  is 
the  incident  field  equivalent  admittance  computed  using  (51).  The  latter  is  due  to  the 
transformed  surface  discontinuity  from  (30).  The  sign  of  Kqq  is  negative  since  it  is  due  to  an 
incoming  wave.  The  same  is  true  for  the  computation  of  Y'qq  in  the  denominator  of  the  first 
term-consequently  it  will  be  the  complex  conjugate  of  the  corresponding  Yqo  in  the  second 
term.  With  those  caveats,  the  equivalent  admittance  for  the  incident  field  arriving  at 
through  a  seeond  layer  is 


y'  y  A/2 

-*  qOO  ^  qOO 


M2 


M2 


/ 


^2^  ^7  ^900  sin(A'( 


M2 

00 


^2) 


(58) 


2.7.  Transmission  and  Reflection  Coefficients 

The  two  types  of  analyses  shown  in  Figure  1,  array  antennas  and  periodic  screens, 
give  rise  to  two  types  of  reflection  terms:  waveguide  mode  and  Floquet  mode.  Waveguide 
reflection  terms  are  computed  in  the  same  way  regardless  of  whether  there  are  superstrate 
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layers.  Those  waveguide  terms  are  discussed  adequately  in  [5]  and  [6].  The  following 
derivations  apply  to  the  free-space  Floquet  mode  transmission  terms  for  both  analysis  types, 
as  well  as  the  reflection  terms  for  periodic  screens. 


2.7.1.  Free  Space 

The  excitation  coefficients  Cp„„  for  Floquet  modes  above  the  structure  (in  the  z  >  0 
halfspace)  are  given  by  (17)  as  the  products  of  the  transverse  field  on  Fg  and  the  mode 
function,  integrated  over  the  unit  cell  cross  section.  These  are  proportional  to  the 
transmission  coefficients:  when  the  incident  and  transmitted  modes  are  different,  the  square 
root  of  their  admittance  ratio  must  be  included.  Let  Tp„„  denote  the  transmission  coefficient 
in  the  +z  direction  for  mode p,m,n: 


+ 

T 

pmn 


\ 


pmn 


er. 


pmn 


(59) 


The  summation  is  an  approximation  to  the  integral  in  (17),  with  the  electric  field 
approximated  using  finite  elements.  For  antenna  radiation,  is  the  admittance  of  a 
waveguide  mode.  For  periodic  screen  reflection  and  transmission  it  is  the  admittance  of  the 
0,0-order  Floquet  mode  in  free  space. 

For  periodic  screens,  the  modal  excitation  coefficients  for  waves  reflected  into  the  -z 
half  space  are  given  by  (20).  The  reflection  coefficients  for  all  except  the  lowest  order 
mode  are  given  by  (59)  with  Fa  replacing  Fg.  The  lowest  order  coefficient,  the  reflectivity 
is 


pOO 


(60) 
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where  q  represents  the  incident  wave  polarization. 


2.7.2  With  Exterior  Dielectric  Layer 

The  reflection  and  transmission  coefficients  must  be  found  from  the  inner  products 
of  the  mode  functions  with  the  transverse  fields  over  the  boundaries  and  Fg  although  the 
desired  response  is  that  measured  outside  the  exterior  dielectrics.  The  modal  excitation 
coefficients  are  the  Cp„„'s  in  (36)  and  (37),  and  it  is  necessary  to  relate  them  to  the  inner 
products  calculated  in  terms  of  Ap^„  and  Bp^„.  If  (38)  is  evaluated  at  z  =  0  and  then  the  inner 
products  of  both  sides  with  each  mode  function  are  taken,  the  result  is 

{a  +B  )  =  f  e' 'glmn  dy 

\  pmn  pmn  /  /  c>pmn  (ol) 

unit  cell 


Since  the  existing  computer  code  calculates  the  inner  product  on  the  right,  all  that  is 
necessary  now  is  to  find  Cp^„  in  terms  of  (Ap„„+Bp„„).  Rewriting  (40): 


C.„=cos(/cjr„0(^ 


+  B  ) +jsm(K^  t)  {A 

pmn  pmn  ^  ^  mn  ^  ^ 


pmn 


B  ) 

pmn  ' 


pq  mO  «0 


(62) 


The  reflected  and  transmitted  terms  only  from  (50)  are: 


./ 


{A  -B  )  = 

^  pmn  pmn^ 


pmn 

j  A/ 
pmn 


(A  +5  ) 

pmn 


(63) 


and  substituting  (29)  and  (63)  into  (62)  leads  to 


C  =  -  6  d  ^6  +  {A  +B  ) 

pmn  pq  w  0  kO  pmn  pmn^ 


pmn 


^pm«cos(ic^„0+yF  sin(x-j^„0 

pmn  mn  '  *'  pmn  ^  mn  ' 


(64) 
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which  is  the  desired  result.  Therefore,  the  reflection  coefficients  for  the  -z  half  space  are 


pmn 


\ 


pmn 


/ 

V 

pmn 


E  c; 


•er. 


h  - S  S  ^3 , 

pmn  pq  mO  «0 


(65) 


where  is  the  admittance  ratio  appearing  on  the  right  side  of  (64): 


y 


/ 

pmn 


M  ^  pmn  .  ^  _M 


cos(/r^„0  +j 


.,M 

pmn 


%  /  iVl  V 


-1 


(66) 


The  incident  mode  admittance,  is  the  free  space  modal  admittance,  as  given  by  (35).  The 

transmission  coefficients  for  the  +z  half  space  are  also  given  by  (59),  without  the  Kronecker 
delta  term  due  to  the  incident  field. 


2.7.3.  Multiple  Dielectric  leavers 


Extension  of  the  above  to  multiple  layers  is  not  a  straightforward  matter  of  using  (66) 
recursively.  Instead,  the  formula  must  be  derived  by  matching  fields  across  a  second 
boundary  in  a  similar  manner  to  Section  2.6.3. 

Let  the  fields  in  the  innermost  of  two  layers  be  represented  as  in  (52)  and  (53).  Then 
matching  fields  across  the  boundary  between  the  two  dielectric  layers  gives  (54),  which  may 
be  substituted  into  (64)  to  give 


pmn 


pq  mO  wO 


■y 


pmn 


a  e 

pmn 


.  M2 


b  e 

pmn 


M2 


(67) 


Expanding  the  complex  exponentials  gives  a  form  similar  to  (62): 
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(68) 


C  =-S  S 

pmn  pq  mO  nO 


'^ypmn 


M2 


M2 


cos(K„L)(a  +6  ) +ysin(/f  „  L)  (a  -b  ) 

^  mn  2^'“  pmn  pmn  '  ''  ^  mn  2-'  ^  pmn  pmn  ^ 


By  analogy  from  (63),  {a  - 

^pmn  )  in  terms  of(«_  +  Uis 


y" 

(  ^pmn  -  ^pmn)  =  "  (  ^pmn  ^  ^ pmn) 

^ pmn 


(69) 


and  substitution  of  this  along  with  (30)  using  the  recursion  substitutions  from  Section  2.5 
into  (68)  gives 


C  =  -  3  3.3 

pmn  pq  m\J  nO 


/ 


+  (d  +  b  )y 

pmn  pmn^  pmn 


M2 

pmn 


^M2  .  M2  ^  ^  ^  ^  \ 

^p«„cos(;c^„r2)+;r  sin(A;„ 


(70) 


where  Y’p„„  is  computed  according  to  (29).  This  is  the  necessary  form,  giving  the 
coefficients  Cp„„  in  terms  of  {ap„„  +  bp„„),  which  take  the  place  of  (Ap„„  +  Bp„„)  in  the  inner 
product  of  (47).  The  quotient  in  (70)  is  a  second  admittance  ratio,  which  may  be  denoted 
y”p„„.  Thus,  the  approach  for  transmission  and  reflection  coefficients  for  multiple  layers  is 
not  to  use  a  formula  recursively,  but  to  multiply  their  admittance  ratios  together. 

2.8.  Summary  of  Extensions  to  Theory 

The  most  general  case  in  terms  of  the  foregoing,  the  one  that  embodies  all  elements 
of  the  extensions  to  theory,  is  a  periodic  screen  that  has  external  dielectrics  on  both  sides. 
To  solve  for  the  properties  of  such  a  structure,  the  following  extensions  to  the  theory  need 
to  be  incorporated  in  the  computer  code: 

A.  In  computing  the  matrix  entries  associated  with  and  1^  ,  the  equivalent 
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admittance  from  (30)  is  used  in  place  of  the  admittance  in  (24).  This  effective  admittance 
represents  the  reflection  coefficient  for  a  Floquet  mode  transformed  back  to  the  radiation 
boundary  through  the  dielectric(s). 

B.  In  computing  the  right  hand  side  terms  for  Fa,  an  equivalent  admittance  from  (51) 
or  (58)  is  used  in  place  of  the  admittance  in  (32).  This  effective  admittance  represents  the 
modal  admittance  of  the  incident  field  transformed  forward  through  the  dielectric(s). 

C.  In  computing  the  power  reflection  and  transmission  from  the  fields  (determined 
by  solving  the  matrix  equation),  the  admittance  ratio  from  (66)  must  be  included.  This 
accounts  for  a  transformation  of  the  Floquet  mode  admittance  back  through  the  dielectric(s) 
to  the  integration  surface  where  the  field  solution  is  available. 

For  array  radiation  analyses,  step  (B)  is  skipped  and  (A)  and  (C)  are  carried  out  for  Fg  only. 
For  array  reflection,  all  three  steps  are  also  required. 

This  completes  the  presentation  of  theoretical  derivations.  The  following  chapter 
presents  test  cases  used  to  verify  the  above  formulas. 


3.  DEMONSTRATION  CASES 

3.1.  Circular  Waveguide  Phased  Array  with  Dielectric  Sheath 

Figure  5(a)  shows  the  geometry  of  an  array  of  circular  waveguides  that  open  into  free 
space  through  a  ground  plane.  There  may  also  be  a  dielectric  layer  just  outside  the 
waveguide  openings.  Figure  5(b)  is  a  view  of  the  finite  element  model,  which  represents  a 
short  section  of  a  single  waveguide  just  below  the  ground  plane.  Waveguide  mode  boundary 
conditions  are  imposed  on  the  -z  side  to  simulate  the  effect  of  an  arbitrarily-long,  match- 
terminated  waveguide  as  discussed  in  [5]  and  [6].  The  incident  field  carried  by  the 
waveguides  is  in  the  TEji  mode,  polarized  in  the  y  direction  (q=2).  The  periodic  radiation 
condition  discussed  in  section  2.3  is  imposed  on  the  +z  side. 

Figure  6  compares  HFEM  calculations  with  multimode  results  from  Amitay  [2:299] 
for  the  active  reflection  coefficient  magnitude  for  scanning  in  the  E  plane  (^  =  90°).  The 
three  solid  curves  are  the  multimode  results  for  slab  thicknesses  of  0,  .1  X,  and  .5  X.  The 
three  sets  of  discrete  points  are  the  corresponding  HFEM  results.  The  fact  that  the  two 
results  agree  so  closely  confirms  the  validity  of  the  expressions  derived  for  the  radiation 


(a) 


(b) 


Figure  5.  Circular  Waveguide  Array:  (a)  Lattice  Geometry; 
(b)  Mesh  of  Short  Waveguide  Section 
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Figure  6.  Comparison  of  Multimode  [2]  and  HFEM  Calculations  for  Active  Reflection 
Coefficient  of  a  Circular  Waveguide  Array  with  Dielectric  Covers 


boundary  matrix,  specifically,  the  effective  admittance  given  by  (30).  This  test  case  also 
validated  the  accuracy  of  the  equivalent  admittance  (66)  used  for  computing  the  transmission 
and  reflection  coefficients  because  the  total  power  transmitted  into  all  propagating  Floquet 
modes  plus  that  reflected  back  into  the  waveguide  summed  to  unity  for  all  angles. 

Figure  7  shows  the  number  of  iterations  required  to  solve  the  matrix  in  each  of  the 
three  cases  fi-om  Figure  6  using  the  biconjugate  gradient  method  (BCGM)  [7].  Consistent 
with  prior  experience,  the  number  of  iterations  increases  near  scan  blindnesses.  Otherwise, 
it  is  nearly  the  same  for  all  three  cases.  The  addition  of  the  exterior  dielectric  has  very  little 
effect  on  convergence,  and  therefore  very  little  effect  on  the  matrix  condition  number.  This 
substantiates  the  claim  made  in  the  introduction  that  the  additional  capability  for  external 
dielectrics  comes  at  almost  no  added  cost  in  solution  time. 
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Figure  7.  Number  of  Bieonjugate  Gradient  Iterations  Needed 
to  Solve  Matrices  for  Cases  in  Figure  6 


3.2.  Rectangular  Waveguide  Array  with  Offset  Dielectric  Layer  (Radome) 

3.2.1.  Radiation 

Figure  8  shows  the  geometry  of  an  array  of  open-ended  rectangular  waveguides  with 
an  air  gap  between  the  ground  plane  and  a  dielectric  layer,  representing  an  antenna  with  a 
radome,  and  an  air  gap  between  the  radome  and  the  radiating  aperture. 

Figure  9  shows  the  grids  for  two  unit  cell  models  to  contrast  the  ways  HFEM  can  be 
used  to  analyze  the  radiation  from  this  array.  In  the  first.  Fig.  9a,  the  mesh  model  represents 
a  section  of  a  unit  cell  outside  the  waveguide  opening,  including  both  the  air  gap  and  the 


28 


DIELECTRIC  WAVEGUIDE 
AIR  GAP  AP^«TURE 


.05  A 


Figure  8.  Geometry  of  Rectangular  Waveguide  Array  with 
Air  Gap  and  Dielectric  Layer:  (a)  Top  View;  (b)  Side  View 


Figure  9.  Mesh  Models  for  Rectangular  Waveguide  Array:  (a)  Unit  Cell 
Section  Above  Aperture;  (b)  Waveguide  Section  Below  Aperture 
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dielectric  layer.  The  shaded  area  is  the  ground  plane’s  conducting  surface  outside  the 
radiating  aperture.  Periodicity  conditions  are  imposed  on  the  unit  cell  side  walls  parallel  to 
z.  That  is  the  approach  that  was  previously  necessary.  In  contrast,  Fig.  9b  shows  the  new 
alternative,  in  which  the  mesh  is  a  single  layer  of  cells  inside  the  waveguide. 

Figure  10  is  a  comparison  of  the  calculated  results  for  active  reflection  coefficient 
magnitude  using  both  models,  with  the  “radome”  layer  having  permittivity  of  1.0  and  4.0. 

In  each  case,  the  incident  field  is  the  TEjo  mode,  whose  electric  field  is  oriented  in  the  y 
direction.  The  fact  that  these  calculations  agree  confirms  the  validity  of  the  recursion 
formula  (30)  for  effective  admittance  used  in  the  matrix  fill  computations. 

Figure  1 1  compares  the  magnitude  of  the  active  element  gain  versus  scan  angle  for 
the  two  cases.  The  agreement  between  these  results  establishes  the  validity  of  (70),  the 
transmission  coefficient  for  multiple  layers.  The  phases  of  both  the  reflection  and 


Figure  10.  Active  Reflection  Coefficient  Magnitude  Calculated  Using 
Alternative  Models  from  Figure  9 
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Figure  1 1.  Active  Element  Gain  Calculated  Using 
Alternative  Models  from  Figure  9 

transmission  coefficients  exhibit  similar  behavior,  although  with  offsets  appropriate  to  the 
difference  in  reference  locations.  The  incident  field  is  referenced  to  Fa,  which  is  different 
in  the  two  cases.  The  transmitted  field  phase  reference  is  at  Fb  (.  1  A,  above  the  ground  plane) 
for  the  grid  in  Fig.  9a.  For  the  grid  in  Fig.  9b,  the  transmission  phase  is  referenced  to  the 
outer  surface  of  the  dielectric  when  one  is  present;  or  to  the  ground  plane  when  it  is  absent. 

The  table  below  compares  the  two  models’  sizes  and  execution  times.  The  new 
method  results  in  a  factor  of  5  time  savings  for  this  problem.  Furthermore,  the  mesh  density 
was  the  same  in  each  case,  appropriate  to  the  required  density  in  the  dielectric  layer.  That 
density  is  excessive  for  the  waveguide  model,  which  contains  free  space  cells  only.  In 
situations  involving  thin,  dense  dielectrics,  the  results  would  even  more  dramatically  favor 
the  new  approach. 
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Table  1.  Comparison  of  Size  and  Execution  Time  for  Two  Models  of  Rectangular 
Waveguide  Array  Antenna  with  Offset  Dielectric  Layer 


Model 

Mesh 
Edges  (N) 

Edges  in 
Port  A 

Edges  in 
Port  B 

Average 

Iterations 

Average 
CPU  secs. 

Unit  Cell  (9a) 

3841 

381 

690 

.112N 

354 

Waveguide  (9b) 

1449 

381 

381 

.116N 

70 

3.2.2.  Reflection 

To  validate  the  formulas  for  the  incident  field  arriving  through  an  exterior  dielectric 
layer,  the  array  reflection  case  is  examined.  The  receive  case  is  modeled,  in  which  the 
incident  field  is  a  plane  wave  impinging  on  the  array  face  fi-om  the  -z  half  space.  By 
reversing  the  orientation  of  the  model  in  Figure  9a,  it  may  be  used  to  simulate  reflection  from 
the  same  waveguide  array  antenna.  The  model  in  Figure  9b  remains  the  same,  with  the 
boundary  conditions  on  each  side  reversed.  The  former  case  is  the  reference  solution  since 
it  uses  the  pre-existing,  validated  methods. 

Figure  12  compares  the  active  element  gain  (power  coupled  into  the  waveguide)  using 
the  two  different  mesh  models,  for  scanning  (direction  of  incidence)  in  both  the  E  (y-z)  and 
H  (x-z)  planes.  As  expected,  the  E  plane  results  are  identical  to  those  in  Figure  11.  The 
correspondence  between  the  two  models  indicates  that  the  recursive  formula  for  the  incident 
field  (58)  is  correct. 

This  concludes  the  discussion  of  computer  code  validation.  The  three  cases  presented 
above,  circular  waveguide  array  radiation,  rectangular  waveguide  array  radiation,  and 
rectangular  waveguide  array  reflection,  verify  the  three  essential  developments  that  were 
needed  to  extend  the  existing  theory  to  handle  all  configurations  of  stratified  exterior 
dielectrics. 
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ACTIVE  ELEMENT  GAIN 


Figure  12.  Received  Power  (Active  Element  Gain)  Calculated  Using 
Alternative  Models  from  Figure  9  for  Plane  Wave  Incidence 
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4.  CONCLUSIONS 


A  method  for  incorporating  the  effects  of  exterior  stratified  dielectric  layers  in  hybrid 
finite  element  calculations  for  periodic  structures  has  been  developed  and  demonstrated. 
This  extension  to  the  method  makes  it  more  efficient  for  modeling  antennas  that  have  wide- 
angle  matching  layers  and/or  radomes.  It  was  previously  necessary  to  extend  the  finite 
element  mesh  through  all  such  layers,  adding  to  the  problem  size  and  execution  time. 

Three  modifications  to  the  existing  theory  were  required.  First,  the  modal  admittances 
used  in  computing  the  matrix  terms  for  a  periodic  radiation  boundary  were  replaced  with 
effective  admittances  that  account  for  the  discontinuities  between  dielectric  layers,  and  the 
layer  thicknesses.  Second,  a  similar  effective  admittance  was  derived  for  the  incident  field, 
which  affects  the  calculation  of  the  right-hand-side  vector  in  the  matrix  equation.  Last,  after 
the  matrix  solution  yields  the  coefficients  for  the  electric  field,  the  transmission  and 
reflection  coefficients  are  computed  using  a  modified  formula  that  includes  multiplication 
of  admittance  ratios  for  each  layer. 

The  foregoing  extensions  to  the  theory  were  demonstrated  to  be  valid  using  an 
updated  version  of  an  existing  hybrid  finite  element  computer  code.  Calculations  of  active 
reflection  coefficient  and  active  element  gain  for  both  transmitting  and  receiving  waveguide 
antenna  arrays  were  shown  to  be  correct.  These  calculations  also  demonstrated  a  factor  of 
five  savings  in  execution  time  for  an  array  antenna  with  a  radome. 
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APPENDIX  A 


NUMERICALLY  STABLE  ADMITTANCE  COMPUTATION 

Calculation  of  the  effective  admittance  for  external  dielectric  layers  involves  terms 
of  the  form: 


Y'  _  cos(2)  +jAmi{z) 
-(4  cos(z) +ysin(z) 


(71) 


where  z  =  x  +jy  and  A  =  Y^/Y^.  When  z  is  large  (in  magnitude)  and  complex,  which  will 
be  true  for  large  orders  of  mode  indices  m  and  n,  the  complex  cosine  and  sine  functions  may 
cause  floating  point  overflow  because  they  involve  hyperbolic  functions.  To  avoid  that 
difficulty,  the  identities: 


cosz  =  cosx  coshj^  -ysinx  sinhj; 
sinz  =  sinx  coshj;  -ycosx  sinhy' 
are  used  to  rewrite  (18)  as 


cos(x) 

l-^tanh(>;)]  +7  sin(x)  -tanh(y) 

cos(x) 

A  -tanh(y)]  -j  sin(x)[l  -^tanh(y) 

(72) 

(73) 


(74) 


In  computing  the  hyperbolic  tangent,  one  of  the  following  two  forms  is  used: 


tanh(>;)  =  ^ — 

+e 


1  -  g 

I  + 


e^y-1 

+  I 


choosing  (22)  forj^O  or  (23)  fory<0. 


(75) 

(76) 
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